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Abstract. The fluid-scale evolution of relativistic magnetic reconnection is investigated by using 
two-fluid and magnetofluid simulation models. Relativistic two-fluid simulations demonstrate the 
meso-scale evolution beyond the kinetic scales, and exhibit quasisteady Petschek-type reconnec- 
tion. Resistive relativistic MHD simulations further show new shock structures in and around the 
downstream magnetic island (plasmoid). General discussions on these models are presented. 
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INTRODUCTION 

Magnetic reconnection is an important fundamental process in plasmas. It has drawn 
growing attention in extreme astrophysical sites such as pulsars and magnetars, where 
plasmas primarily consist of electron-positron pairs. One of the most notable applica- 
tions is the magnetic dissipation problem in a pulsar wind, an ultrarelativistic plasma 
flow from the pulsar magneto sphere. It is expected that an equatorial current sheet (like 
the heliospheric current sheet) flaps extremely due to the oblique rotation of the central 
neutron star, and that reconnection inside the current sheet dissipates the magnetic en- 
ergy [1]. Importantly, in those environments the magnetic energy exceeds the rest mass 
energy of lightweight electrons and positrons. When reconnection transfers magnetic 
energy to that of plasmas, we need to take into account special relativity effects both in 
the bulk motion and in the plasma heat. 

Our understanding of relativistic reconnection is much more limited than our under- 
standing of the nonrelativistic counterpart. In the last decade, there has been some fun- 
damental work in the following two areas: magnetohydrodynamic (MHD) theories of a 
steady-state structure [2, 3, 4] and time-dependent particle-in-cell simulations on kinetic 
scales [5, 6, 7, 8, 9]. Since these two are temporally and spatially separated, it has been 
difficult to associate the results from the two research areas. After an initial attempt [10], 
MHD simulations of the basic reconnection process have been stalled for a while. 

In order to bridge the two areas, we have developed two simulation models to study 
relativistic magnetic reconnection: the electron-positron two-fluid model [11, 12] and 
the resistive relativistic MHD (RRMHD) model [13]. 



TWO-FLUID SIMULATIONS 



The basic equations consist of relativistic fluid equations for electrons and positrons, and 
Maxwell equations [11, 12]. For simplicity, c is set to 1. 
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In these equations, the subscript p means positron properties (and e for electrons), 7 
is the Lorentz factor, n is the proper density, m = yv is the 4-vector, w is the enthalpy, 
p is the proper pressure, and qp — —qe is the charge. We assume a F-law equation 
of state with the adiabatic index F = 4/3. Therefore the enthalpy is given by w = 
nmc^ + [F/ (F — \ )]p. In order to mimic the effective resistivity, we introduce an inter- 
species friction force in the last terms of the momentum and energy equations with a 
coefficient Tfr. Internally, we consider the energy density relative to the rest mass energy 
density (i.e., we subtract Eq. (1) from Eq. (3)) similarly as in earlier work [10, 14]. 

The time evolution is solved by a standard numerical scheme, i.e. a modified Lax- 
Wendroff scheme. A common difficulty of relativistic fluid simulations is that the fluid 
macro properties (the "conservative variables") on the left hand sides are nonlinear 
combinations of the basic elements (the "primitive variables") such as 7 and p [14]. In 
our simulation, we calculate the primitive variables from the conservative variables by 
analytically solving a quartic equation of |m| [11] on all grid cells at each half timestep. 

We investigate a two-dimensional system evolution in the x-z plane. The reconnection 
point is set to the origin (x, z) = (0, 0) . We employ a Harris-like initial configuration: B = 
Bq tanh(z) jc, j = Eq cosh~'^(z) y,n = no cosh~^(z) -|- 0. Ino^ P = nmc^, « = 0, and £ = 77 j. 
In the upstream region (|z|3>0), the magnetization parameter is a = BQ/[4n{2w)] = 4. 

The relevant Alfven speed is CA,up = [(7/(1 + c)]^/^ ~ 0.89c or YCA,up = = 2. We 
assume that the effective resistivity t] oc Zfr is localized around the reconnection point. 
Neumann-like boundaries are located at x = ±120 and at z = ±60. 

Reconnection takes place around the center. Figure la shows the x-component of the 
4-velocity, Ux = Y^x, at ? = 100 in unit of the light transit time c^^. Magnetic field lines 
are transported from the background regions, cut and reconnected at the reconnection 
point, and then ejected with the reconnection jets. One can see the bi-directional jets 
from the reconnection point and magnetic islands (plasmoids) in front of the jets. The 4- 
velocity of the reconnection jet is typically Ux ~ 2, which is Alfvenic with respect to the 
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FIGURE 1. (Color online) The plasma 4-velocity Ux and the magnetic field lines at (a) t ~ 100 and (b) 
f = 400. In Panel (a) the arrows indicate the vertical structures, which will be discussed later. In Panel (b) 
the dashed fines represent Petscfiek-type slow-sfiock regions. (From run L3 in Ref. [11]) 



upstream condition, and it is ~ 3.5 at the local maximum right behind the plasmoid. 
Reconnection sets up an inflow speed of ~ 0.14-17. Since the reconnection rate or a 
normalized form of the flux transfer speed is an order of ^ ~ iyinlvoui) ~ 0(0. 1), this 
is a fast reconnection. Later, the plasmoids reach and go through the open boundaries at 
X = ±120. Although there are minor reflections from the boundaries, the system exhibits 
a quasisteady structure in the long term. Figure lb presents the late-time snapshot at 
t = 400. The reconnection jets are confined by slow-shock-like regions, as indicated by 
the dashed lines in Fig. lb. One can see a quasisteady Petschek-type structure [15]. 

Interestingly, our analysis reveals that the Petschek outflow becomes narrower than 
the nonrelativistic counterpart, as theoretically predicted by Lyubarsky [4]. Meanwhile, 
reconnection proceeds faster, especially in the ultrarelativistic, high-a, regimes. It was 
also found that an out-of-plane magnetic field (guide field) drastically changes the 
composition of the energy outflow. Without a guide field the main energy carrier is the 
plasma enthalpy flux, while the Poynting flux carries the energy in the presence of a 
moderate guide-field [12]. 



RESISTIVE MHD SIMULATIONS 



We use the following RRMHD equations [10, 16] in Lorentz-Heaviside units with c = 1. 

^r(rp)+v-(pM) = 0, (8) 

dt{rwu + ExB)+V-i^{p^ — )I + wuu-BB-EEj = 0, (9) 

dt{fw-p^ ^ ) + y-{rwu + ExB) = 0, (10) 

dtB + VxE^O, dtE-VxB = -j, dtPc + y-j = 0, (11) 
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Here, pc is the charge density. We developed three RRMHD codes [10, 16, 17]. Among 
them, the present simulations are carried out by the HLL-type code [16]. The RRMHD 
schemes differ from nonrelativistic schemes in the following two ways. One is related to 
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FIGURE 2. (Color online) (a) The x-component of 4-velocity (m^) at the developed stage. The contour 
lines indicate the magnetic field lines, (b) The m, -profile in the front side of the plasmoid. The domain is 
identical to the yellow box in Panel (a). (From run 1 in Ref. [13]) 



the relativistic primitive variables. We use the same quartic equation solver to recover the 
primitive variables. The other involves the usage of Ampere's law. In the nonrelativistic 
MHD, we assume 7 = V x 5 and then derive E from Ohm's law. In the relativistic 
MHD, we use the Ampere's law to advance E. Configurations are similar to those in the 
two-fluid case. We employ left-right symmetry to reduce computational cost. We use a 
localized resistivity rj = t](x,z) in the relativistic Ohm's law (Eq. 12). 

The system evolves similarly in the two-fluid and MHD runs. In the present case, 
we can resolve sharper shock structures because we employ a shock-capturing scheme. 
In the two-fluid model, the system contains the Larmor radius or the inertial length 
and so it makes no sense to discuss shorter structures. We have no such concern in the 
scale-free RRMHD model. Figure 2a shows the Wx -profile in the developed stage. This 
is a Petschek-type reconnection geometry. The typical outflow is Alfvenic, ~ 2. The 
jet hits the plasmoid at a fast shock ("FS" in Fig. 2a) at the maximum 4-velocity of 
Ux ~ 2.75. 

The simulation reveals new shock structures: postplasmoid vertical slow shocks, for- 
ward vertical slow shocks, and the shock-reflection structure. One can see the postplas- 
moid shocks at X ~ 60, as indicated by "SS-1" in Fig. 2a. The shock propagates in the 
-|-x-direction at the front of the reverse plasma flow {Ux < 0; blue regions in Fig. 2a). 
As the plasmoid moves, it compresses the surrounding plasmas on both upper and lower 
sides, and the cavity region also appears behind the plasmoid. The pressure gradient be- 
tween those two regions drives the reverse flow along the field lines [13, 18]. The right 
side is the shock upstream. One can also recognize similar shock-like structure in the 
two-fluid run, as indicated by arrows in Fig. la. We find another small shock outside the 
plasmoid ("SS-2" in Fig. 2a; see also Fig. 2c in Ref. [13]). This is a relativistic version 
of the forward vertical slow-shock [18]. We further analyze the properties of these two 
shocks. In the shock normal direction [19], we evaluate the shock speed v shock and the 
MHD wave speeds in the lab-frame [20]. Key results are presented in Table 1. In both 
cases, on the right (upstream) sides, the shocks travel faster than the outgoing slow- 
mode, but slower than the Alfven mode, Vs+,r < Vshock < ca+,r- This is the character of 
slow shocks. Finally, we find the shock-reflection structure in the front side, as indicated 
by the small box in Fig. 2a. Figure 2b shows the Wx-profile in the same domain. This 
one takes advantage of the current sheet configuration. Since there is initially a dense 



TABLE 1. Selected shock properties. 





Location 




^shock 


Vs+,L 




CA+,R 


SS-1 
SS-2 


(6L5, 5.0) 
(105.9, 8.2) 


(1.00,-0.08) 
(0.90,0.43) 


0.41 
0.62 


0.50 
0.66 


0.29 
0.56 


0.80 
0.80 



The shock normal vector h, the shock velocity Vshoch the slow-mode, and the Alfven wave speeds on the 
left (L) and right (R) sides of the shock. All velocities are evaluated in the lab-frame in the n-direction. 



plasma near the center z ~ 0, plasma flows are bifurcated. The twin flows and the central 

immobile plasmas are separated by intermediate shocks [21] (the dashed lines or "IS"). 
As seen in Fig. 2b, the twin edges invoke oblique bow shocks (fast shocks: "FS" in Fig. 
2b), and then they are reflected by the intermediate shocks many times. In the long term, 
the system exhibits a chain of diamond-shaped structures. Inside the central immobile 
region, the magnetic fields are weak or zero. The fastest wave is the sound speed in 
the current sheet Cs,cs ~ 0.53. The twin edges move faster ~ 0.75 > Cs,cs, because they 
are essentially driven by the reconnection jet ~ CA,up ~ 0.89. We therefore expect the 
obUque shocks when CA,up > Cs,cs- 

DISCUSSION 

These simulations have revealed new insights in the fluid-scale properties of relativistic 
reconnection. The relativistic Petschek reconnection is fast, it features an Alfvenic out- 
flow jet, and the outflow exhaust becomes narrower and narrower as the magnetization 
parameter a increases. The overall results confirm Lyubarsky's theory [4]. In our simu- 
lations, (7 ranges up to ~10. Extended parameter surveys and many basic issues need to 
be investigated. 

The simulations exhibit new shock structures, especially in the RRMHD model. We 
think that they appear when the Alfvenic reconnection jet exceeds the sound speeds 
[13, 18]. Comparing the outflow speed ^ CA,up = [07/(1 + c)]^/^ and the upper limit of 
the sound speed 1 / \/3, a sufficient condition for a supersonic reconnection jet is 

o > 1/2. (13) 

Therefore we expect various MHD shocks in the reconnection system in high-c regimes. 

At present, the most important issue is the effective resistivity, which crucially con- 
trols the system evolution. In this work, we have employed a spatially-localized resis- 
tivity for our main results. It is known that such resistivity leads to a Petschek-type fast 
reconnection in the nonrelativistic MHD studies [22]. Meanwhile, other resistivity mod- 
els give significantly different pictures: a turbulent outflow with secondary islands [11], 
a slow Sweet-Parker reconnection and so on [13]. We need to find a practical resistiv- 
ity model to reproduce realistic reconnection, by referring to the kinetic results [7]. Of 
course a true resistivity is desirable, but this is a long-standing problem in reconnection 
physics. 

The RRMHD model has another problem to overcome. The standard equations are 
very stiff when the magnetic Reynolds number 5 ~ 1 /t] is high. Implicit-type schemes 



[17, 23] have been recently developed to explore high-i' regimes and applications to 
Sweet-Parker reconnection are in progress [24]. 

The development of the two-fluid and RRMHD models will benefit the large-scale 
modeling of relativistic MHD instabilities [25, 26]. Even in ideal MHD studies, the 
late-time evolution often exhibits potential reconnection sites with antiparallel magnetic 
fields or (numerical) reconnection. Reconnection with a physical resistivity will be a 
basic element of these important problems. 

Our results have implications for space physics also. Inspired by these results, we 
carried out nonrelativistic MHD simulations and found the same shock structures in 
low-beta plasmas [18]. In other words, our discoveries are ubiquitous in both relativistic 
and nonrelativistic reconnections. In the solar corona, the fast-shock-reflection structure 
at the plasmoid front may be related to energetic particle acceleration. 
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